Draft version February 2, 2008 

Preprint typeset using LATpX style emulateapj v. 04/03/99 



O 
O 

(N 

u 
< 



O 
H 



1. INTRODUCTION 



WIDE ANGLE RED SHIFT DISTORTIONS REVISITED 

ISTVAN SZAPUDl" 

"Institute for Astronomy, University of Hawaii, 2680 Woodlawn Dr, Honolulu, HI 96822 

Draft version February 2, 2008 

ABSTRACT 

We explore linear redshift distortions in wide angle surveys from the point of view of symmetries. We show 
that the redshift space two-point correlation function can be expanded into tripolar spherical harmonics of zero 
total angular momentum <S , i 1 / 2 ; 3 (£i,X2, x). The coefficients of the expansion Bj 1 i 2 ; 3 are analogous to the C;'s of 
the angular power spectrum, and express the anisotropy of the redshift space correlation function. Moreover, only 
a handful of S/ 1 ; 2 ; 3 are non-zero: the resulting formulae reveal a hidden simplicity comparable to distant observer 
limit. The B; 1 i 2 j 3 depend on spherical B esse l moments of the power spectrum and / = n°- 6 /b. In the plane 
parallel limit, the results of Kaiser ( 119871) and lHamiltonl dl993l) are recovered. The general formalism is used to 
derive useful new expressions. We present a particularly simple trigonometric polynomial expansion, which is 
arguably the most compact expression of wide angle redshift distortions. These formulae are suitable to inversion 
due to the orthogonality of the basis functions. An altern ative Legendre poly nomial expansion was obtained as 
fN| ' well. This can be shown to be equivalent to the results of Szal av et"all ( fl9 98). The simplicity of the underlying 

theory will admit similar calculations for higher order statistics as well. 

Subject headings: cosmic microwave background — cosmology: theory — methods: statistical 

It has been known for decades that the two-point correlation function, or power spectrum, measured in redshi ft surveys is d istorted 
■"nT" ' by the peculia r velocities of galaxie s. The anisotr opy of the correlation function was demonstrated by iDavis & Peebles! dl983l) : 
, iPeeblesI dl980l) . In a seminal work, Kaiser (1987) demonstrated that in the plane parallel limit there is a simple transformation 
■ between the redshift space and real space density contrast. This results in an anisotropic enhancement of the power spectrum by 
(l+//i 2 ) 2 , where /i is the cosine of the angle between the line of sight and the wave-vector. This simple formula has become the start- 
C"| . ing p oint of many extensions, which have used expansion into Legendre polynomials (e.g.. lHamiltonin"993tlHam ilton & Culhanel 
Q^' 119961) . or numerical methods (Zaroubi & Hoff man! 1 19961) . Most analyses assume a small opening angle dCole et all Il995t) . i.e. 
they stay essentially in th e distant observer limit. Others works used a expansion with formally infinite number of coefficients 
( iHeavens & Tayloil|1995l). Num erous galaxy re dshift surveys have be en suc cessfully analyz ed with such methods, most nota bly th e 
PSCz lTadros etaTTll999l> . 2dF jPeacock et all PiOOll iHawkins et all l2Q03t iTegmark etliil Eool and SDSS JZehavi et allEool . 
, For a review of methods in the above spirit and the corresponding applications. see lH^rrm7onT(ll998l) . 

To address the needs of wide angle redshift surveys, full treatment of wide angle distortions have been given by Szala v et"af] 
( 1998), where they identify the coordinates in which the expression of the redshift space two-point correlation function is compact, 
and most importantly finite. They have also argued, that the power spectrum will necessarily have an infinite expansion, as it arises 
' from the convolution of the density field with a non-compact kernel. They concluded that correlation functions are more convenient 
for redshift space analyses then power spectra. Their results is suitable and has been used for "forward" analyses, such as the 
Karhunen-Loeve method, in which the cor relation fu nction is predicted and c ontrasted with data. Applications to the SDSS are 
presented most recently by P ope et all ( 120041) . (see also lTegmark et alll2003albl) . 

In this paper we analyze the symmetries of redshift space distortions. The next section shows that zero angular momentum tripolar 
functions form a natural basis to expand the redshift space correlation function, and that only a surprisingly small set of expansion 
coefficients are will be non-zero. In section 3 we present the most important properties of the basis functions, and the connection 
with the Kaiser-Hamilton limit. Section 4 employs the general theory to obtain compact expressions for the redshift distortions using 
conveniently chosen variables. In the final section we present discussions, and conclusions. 

2. REDSHIFT DISTORTIONS IN LINEAR THEORY 

The theory of redshift distortions is based on the redshift to real space transformation, Si = Xi — fvjiiij, where Xi is a unit vector 
pointing to the galaxy from the origin, f = f2°' 6 /6, and Vj is the peculiar velocity in units that its divergence equal to the density 
to linear order (Scoccimarro et allll999l) . The transformation of the density S then can be estimated via the linear theory Jacobian 

J = 1 — fxiXjdjVi, as 

S s (x) = J ^Le ik »*>(l + f(x s k s f)8{k) (1) 

where the Einstein convention was followed for summation of multiple indices. The effects due to the gradient of the selection 
function and the "rocket effect' ' (lKaiserlfT987i) were neglected. As a consequence, the redshift distorted 2-point correlation function 
£ s (xi, X2) reads 



£ s (x 1 ,x 2 ,x,x) 



= (l + i/) 2 £(z) 



+ |/ 2 / Ake ik£x P 2 (kx 1 )P 2 (kx 2 )P(k), (2) 



+(§/ + If 2 ) I ^e ikix P 2 (kx 1 )P(k) + (1 «-> 2) 

1 f 2 f rf 3 i 
"9 J J (2tt 

where P 2 is the second order Legendre polynomial, P(k) is the linear power spectrum. We have chosen to characterize the correlation 
with the two directions from the line of sight, and introduced x = x\ — x 2 , and the c orresponding unit ve ctor. By construction, these 
three unit vectors are in the same plane. The above formula is in full agreement with Szalay et all ill 9981) . 

A few simple observations are in order with respect to the above function i) it is rotationally invariant ii) £ s (xi,x 2 , —x,x) = 
£ s *(xi,x 2 ,x,x) = £ s (xi,x 27 x,x), since the correlation function is real, iii) £ s (£i, x 2 , x, x) — £ s (x 2 ,xi,x,x) the invariance 
of the correlation function under permutation, iv) the three vectors are constrained to be in the same plane, v) the unit vector x 
is constrained to be between x\ and — x 2 . vi) We can extend the function for x vectors outside this range, with the definition 
£ a (xi, x 2 , x, x) = £ s (xi, —x 2 ,x, x). This leaves Equation[2]formally valid, since Pi{— fx) = (—l) l Pi([i), and the equation contains 
only even Legendre polynomials. 

For a function depending on three directions, the tripolar spherical harmonics expansion (e.g., IVarshalovich et all [7988) is the 
most natural 

Or x) = £ B l l l ^{x){Y h ® (Y l2 ® Y h ) l23 ) LM 

= ^ B MM 23 {x)D L mmi {Yh ® (Y l2 ® Y h ) h3 ) LM , , (3) 

where D\ lM , is the Wigner rotation matrix, and the second line follows from rotational invariance. As a consequence only L = 0, 
i.e. zero total angular momentum, can appear in the above expansion. This means that quite generally, the expansion will only 
contain the functions (Y^ ® (Yi 2 ® Yz 3 )z 2 3)oo with ^23 = h- In order to simplify the final expressions, we will use the following 
functions, which are proportional to the zero angular momentum tripolar functions 

S hh i(xi,x 2 ,x)= ( TO \ m 2 l m j c ii mi {xi)Ci 2m2 {x 2 )Ci m {x) (4) 

m\ ,m2,m 

where ( ^ ^ 2 ' ] is the Wigner 3j symbol, and Ci m — ^4tt/21 + lY/ m are proportional to the spherical functions. 

Expanding Equation|2] and using Pi(xix 2 ) = J2 m Qm(xi)Cf m (x 2 ), e lkx = J2i m ( 21 + l ) il 3l{ kx ) C im{k)C* m (x), and the Gaunt 
integral, one find finds that 

B lM h = B i M = (2/ + 1)6(x) ., ^ l Q i b J ^ FhFh: (5) 

where we have introduced ^j(x) = / ^k 2 ji(xk)P(k), the moments of the power spectrum with spherical Bessel functions, and 
F = 1 + 1 /3/, F 2 = 2 /3/, otherwise 0. From property (iii) it follows that B hhh = B hllh 

It is worth to write this result explicitly, since only a few terms are non-zero due to the initial expression and constraints from 
group theory: 

b ooo (x) =(i + im (x) 

B 022 {x)=B 202 {x) = _ {2f+2f2) ^ 2{x) 
B 22 2[x) = ^p Ux) 

B 22 \x) = $£f 2 U*)- (6) 

These functions, not unlike the C; for the angular power spectrum, form a natural basis for maximum likelihood estimation. They 
can be used as an intermediate step for estimation of cosmological parameters in the linear regime. 

3. PROPERTIES OF THE Sl 1 l 2 l 3 FUNCTIONS 

The Sij 2 i 3 functions are a subset of tripolar spherical harmonics. They form an orthogonal complete basis for expanding spher- 
ically symmetric functions depending on three unit vectors. It should be emphasized, that orthogonality is true only when the unit 
vectors are integrated over the full sphere dflidfl 2 dil3 unrestricted. From the definition it follows that 

(4tt) 3 / 2 

dn 1 dn 2 dn 3 S h i 2h (x 1 ,x 2 ,x)S jlhh (x 1 ,x 2 ,x) = S hjl Si 2j2 S hh === ===■ (7) 

y {M\ + l)(M 2 + 1J(2<3 + J-J 

The plane parallel limit (Kaiseri ri987l) can be obtained by assuming that the first two unit vectors are parallel. Using the properties 
of the Wigner coefficients and spherical functions, it is easy to show that 

Si lh i 3 (xi,xi,x) = ( q q q ]Pi 3 (xi-x). (8) 
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As a special case one can see that independently of the unit vector xi 

Sooo = 1 

Soii{xi,x 2 ) = Sioi(xi,x 2 ) = S m (x ll x 2 ) = ^ 2 i +1 p i( x \ ■ x i) 

(9) 

and all other functions with any zero index are zero. 

It follows from the above properties by simple algebra that in the plane parallel limit X\ ~ x 2 

£'(ii,ai,£,aO = (l + |/ + i/ 2 Ko+ 

-(|/ + |/ 2 )6P 2 (i,3) 

+^/ 2 £ 4 P 4 (l ) 3), (10) 
in agreement with lKaiserl ( ll98l) : lHamiltonl(ll993l) : lHamilton & Culhand i ll 9961) . 

4. COORDINATE SYSTEMS 

Since the three unit vectors are constrained in the same plane (property (iv), the redshift space correlation function depends only on 
the shape and size of a triangle. Consequently, the angular dependence can be parametrized by two angles. Since the S\ x m functions 
are rotationally invariant, we can fix the plane of the vectors, and even rotate one of them to a fixed position in order to obtain useful 
expressions. We explored the following choices A) the z axis is perpendicular to the plane of the vectors, and we fix <f> — for x, 
B) z axis coincides with the third unit vector, z = x, C) it points in between the first two unit vectors z oc x\ + x 2 . For the latter 
choices, we can assume that all the vectors are in the <f> = plane. 

For choice A), S/ 1 j 2 /(xi, x 2 , x] = S'; 1 ; 2 /(7r/2, </>i, ir/2, <fi 2 , ir/2, 0), and one obtains an expansion of the form £ s (</>i , 02 , = 
Em n 2 =o 2 a "i"2 cos(rti0i) cos(n 2 (f> 2 ) + b nin . 2 sin(7ii0i) sm(n 2 <f) 2 ), where the only non-zero coefficients are: 

«oo = (l + 2 -i + 3jf ) £o(x) - (I + 2 -£) &x) + ffi &(x) 
ao2 = »20 = (ir _ &( x ) + €*( x ) 

a 22 = Tg£o{x) - It + ^rw&ix) 

b 22 = h €o( x ) - 21 &( x ) - &( x ) 

(11) 

According to property (v) there is a restriction that </>i < <f> 2 . which, however, can be lifted by symmetry property (iii). For 
a fixed x, the two angles can span the full range of the integration, if there are no restrictions from incomplete sky cover- 
age. Then, the above becomes (double) orthogonal expansion, where the coefficients can be obtained simply by integration, e.g. 

a nin2 oc f^ 2 d(pid(f> 2 ^s( ( f > i, 02: x) cos(2n 1 1 ) cos(2n 2 (f> 2 ) . If the correlation function is measured and binned according to the 
above expression, the four independent coefficients can be obtained either by numerical integration, or by fit; the latter would be 
probably preferable for realistic surveys with incomplete sky coverage. 

The above expression can be inverted easily with a computer algebra package: the variables /, £q> £4 can be expressed analyti- 
cally as a function of the four coefficients. Unfortunately, the analytical expression is too complicated to list here (it is a solution of a 
fourth order polynomial in /), but if needed it can be easily obtained with any computer algebra package, such as Mathematica. 
In practical applications, however, numerical inversion is expected to be more robust. 

The plane parallel limit is cj>i — (f) 2 : reexpressing t he trigono metric functions in Legendre polynomials indeed yields, after 
somewhat tedious calculation, the familiar expression of Hamilt on! dl993l) . 

In coordinate system B) Si 1 i 2 i(xi,x 2 , x) = S'/ 1 ; 2 /(6'i, 0, 9 2 , 0, 0, 0), but the meaning of the two angles is the same as for coordinate 
system A). It can be shown with simple but tedious calculat ion that it re duces to the same expression as above. 

Finally, coordinate system C) is identical to that of Szalav etaU l ll998l) : As we show next, it produces a double Legendre expansion. 
For this choice we have Si 1 i 2 i(xi,x 2 ,x) — S , ; 1 ; 2 /(6', 0,6, it, j , 0) (0 = ir ensures that the z axis is between the first two unit vectors). 
The notation of the variables emphasizes the connection with Szalav "etail d!998h . 

Using the fact that Ci m (9, tt) = (— l) m Ci m {0, 0), inserting the Clebsch-Gordan expansion 

C hmi (x)C hm2 (x) = ( ) ( mi m 2 m 3 ) C hm 3 (x), (12) 

and using property (iii), one can show that in 7713 and m are both even. This means that in that case the S; 1 / 2 / functions have only 
even associated Legendre functions, which can be reexpressed into a finite set of Legendre polynomials. Since there is only a few 
non-zero coefficients, we show the explicit result instead of the general derivation, which is not very illuminating. Finally, the redshift 
space two-point correlation function can be expanded into multiples of two Legendre polynomials, £ s (cos 9 = \x\ , cos 7 = fi 2 , x) = 
Szi ; 2 =o 2 c hhPh {^i)Ph (^2), where the only non-zero coefficients are 

coo = (1 + ¥ + 3f ) &(*) + &(*) + 3 Aux) 
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Iff) 6(X) + ^4(X) 

^(z) 



(13) 



In this form the plane parallel limit is /ij = 1, i.e. Pi(l) = 1. It is simple (although a bit tedious) matter to show that this agai n 
returns the right answer. Again, simple, tedious calculation shows that our expression reproduces the results of Szal av et all fl 998). 
if in their Eq. 15 the typographical error 4/15 — > 8/15 is corrected. 

The above is formally an orthogonal expansion. Due to property (v), 9 < 7 < it — 9 must be satisfied for for any given x. 
According to property (v), however, the range can be extended, and the orthogonality of the Legendre polynomials Pi(fi\) is ensured. 
The coefficients in the above expansion can be obtained by Gauss-Legendre integration of the correlation function. Measuring the 
correlation function in these coordinates would produce a method which would be the closest generalization of the original Kaiser- 
Hamilton method. 

The above form perhaps provides the most natural connection to the plane parallel limit, therefore it can be used to quantify 
deviations from it. Figure 3. plots wide angle to plane parallel ratio of the two most common estimators for redshift distortions: 
the ratio of the redshift space and real space correlation functions, and Hamilton's Q(s) — £2(s)/(3/s 3 J £ (y)y 2 dy — £o(s)), the 
modified quadrupole to monopole ratio. According the Figure, a simple restriction of the opening angles at 9 <, 15 — 20 degrees 
would ensure the accuracy of traditional measurements assuming the plan-parallel approximation. 
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FIG. 1. — The two most common estimators for redshift distortions are compared to their wide angle analogue: §n(s) /£M (thick line), and Q(s) = 
?2(s)/(3/s 3 J £,o{y)y 2 dy — £o(s)) (thin line). The ratio of wide angle to plane parallel prediction is plotted as a function of the half opening angle 9. The 
three sets of lines correspond to three slopes of the correlation function 1.5, 1.75, 2.0, increasing and decreasing for the three thick and thin lines, respectively. In a 
realistic measurement, the final result would be a weighted average over a set of opening angles, represented by the above curves. 



5. CONCLUSIONS AND DISCUSSION 



We have analyzed the underlying symmetries of r edshif t space distortions. As a result, we presented three novel expres sions for the 
redshi ft space two-point correlation functions, Eqs. 1611 ll andll 31 The last of these equations turns out to be identical to Szala vetall 
(1998). 

First, we have shown that the two-point correlation function can be expanded into tripolar harmonics and that rotational symmetry 
restricts the non-zero components to those with L = 0. We have found by direct calculation that the quadratic nature of linear 
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redshift distortions restricts the expansion to five unique functions -B; 1 ; 2 ; of Equation[6] each depending on the distance between the 
two points only. These are analogous to the C;'s of the angular power spectrum, constitute a natural basis for maximum likelihood 
analysis of redshift space data. The ful l machinery of maximum likelihood can be adapted naturally (e.g. Vog elev & SzalavUl996t 
Teg mark et all [19981 iBond et alll2000t) . Details are left for subsequent research. 

The rotational invariance of the basis functions, and the fact that the three unit vectors are constrained into a plane, allows us to 
fix a plane and an orientation within that plane. Using this freedom, we identified two convenient coordinate systems, which corre- 
spond to particular choices of variables. The resulting expressions, Eas.ll llandl3lare especially convenient, since they respectively 
correspond trigonometric and Legendre polynomial expansions in the two angular variables. The first expression is possibly most 
c ompact formula for linear redshift space distortions of the two-point correlation function, the second can be shown to be identical 
to lSzalav et alJ(ll998l) . The orthogonality in these expansions presents an opportunity for applications similar in spirit to the original 
Kaiser-Hamilton method, but fully correct for wide angles. 

This paper deals with the theory of redshift distortions, and thus lays the groundwork for possible future applications. We have not 
discussed practical issues, such as incomplete sky coverage, noise. Clearly, even if these issues are important, a maximum likelihood 
technique to find the parameters of the expansions would be still optimal. On the other hand, incomplete sky coverage and noise will 
cause leakage, pos sible emergenc e of higher / anisotropies. We conjecture that such difficulties for direct methods could be solved 
similarly to Szapu dTet al.l (1200 ll) . where the analogous problem for the angular power spectrum was tackled. Details of practical 
applications are left for subsequent research. 

Our calculation could be simply generalized for the effects of the gradient of the selection function. However, Eauationll3lcould 
not be generalized due to odd associated Legendre functions entering in the expansion. These cannot be re-expressed as Legendre 
p olynomials, thus the final results would not admit a Legendre expansion. This can be explicitly demonstrated from the final results 
of Sz alav et"afl(ll998l) . There appear to be no analogous problems when generalizing Equation^J we conjecture that trigonometric 
polynomial expansion will be still possible, with terms of odd orders appearing. The simplicity of the present theory opens up the 
possibility of generalizations for higher orders, such as three-point correlation functions and cumulant correlators. These calculations 
will be presented elsewhere. 

It is a pleasure to thank Alex Szalay for stimulating discussions. The author was supported by NASA through AISR NAG5-1 1996, 
and ATP NASA NAG5-12101 as well as by NSF grants AST02-06243 and ITR 1 120201-128440. 
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